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Michael Mc Gettrick 
De Brun Centre for Computational Algebra, 
School of Mathematics, National University of Ireland, 
University Road, Galway, Ireland 

Jaroslaw Adam Miszczak 
Institute of Theoretical and Applied Informatics, 
Polish Academy of Sciences, 
Baltycka 5, 44-100 Gliwice, Poland 
(Dated: 14/01/2013 (v. 0.69)) 

We study the model of quantum walks on cycles enriched by the addition of 1-step memory. 
We provide a formula for the probability distribution and the time-averaged limiting probability 
distribution of the introduced quantum walk. Using the obtained results, we discuss the properties 
of the introduced model and the difference in comparison to the memoryless model. 
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I. INTRODUCTION 

During the last few years a considerable research effort 
has been made to develop new algorithms based on the 
rules of quantum mechanics. Among the methods used to 
achieve this goal, quantum walks, a quantum counterpart 
of random walks, provide one of the most promising and 
successful approaches. 

Classical random walks can be applied to solve many 
computational problems. They are used, for example, to 
hnd spanning trees and shortest paths in graphs, to find 
the convex hull of a set of points or to provide a sampling- 
based volume estimation [1 . Today a huge research effort 
is devoted to applying random walks in different areas of 
science. Classical random walks find their application in 
a plethora of areas. This has motivated big interest in us- 
ing a similar model for developing algorithms which could 
harness the possibilities offered by quantum machines. 

Quantum walks are counterparts of classical random 
walks governed by the rules of quantum mechanics [T]- 
|3] and provide a promising method for developing new 
quantum algorithms. Among the applications of quan- 
tum walks one can point out: solving the element dis- 
tinctness [1] and subset finding 5 problems, spatial 
search 6 , triangle finding [7\ and verifying matrix prod- 
ucts [8]. The survey of quantum algorithms based on 
quantum walks is presented in [5]. 

In this paper we introduce and study the model of 
quantum walks on cycle [TUl HI] enriched by the addition 
of 1-step memory [12] . Our main contribution is the cal- 
culation of the probability of finding the particle at each 
position after given number of steps, and the limiting 
probability distribution for the introduced model. We 
also point out the differences between quantum walks on 
cycles with and without memory. 

This paper is organized as follows. In Section [TI] we in- 
troduce the model of a quantum walk on cycle with one- 
step memory. In Section |III| we analyze the introduced 



model using Fourier transform method and we discuss the 
behavior of the time-averaged limiting probability distri- 
bution of the discussed model. Finally, in Section [TV] we 
summarize the obtained results and provide some con- 
cluding remarks. 



II. THE MODEL 

In the model discussed in pj] the space used by a quan- 
tum walk is composed of two parts - 1-qubit coin and 
d-dimensional state space, i.e. H = Hi ® Hd- The shift 
operator in this case is defined as 

1 d-l 

S = ^2^2\c)(c\ ® \v + 2c~ 1 mod d)(v\. (1) 

c=0 t'=0 

Here we adopt this model and extend it with an ad- 
ditional register, referred to as memory register, which 
stores the history of a walk. 

For a quantum walk with one step memory one needs 
a single qubit to store the history. In this case we use a 
Hilbert space of the form H = Hi^Hi^T-Ld, respectively 
for a coin, memory and a position. 

As in the case of a memoryless walk, the first register is 
the coin register and the third register is used to encode 
the position of the particle. The second register stores the 
history of the walk. The history is encoded as direction 
from which the particle was moved in the previous move. 
If this register is in the state |0), the previous position 
of the particle was n + 1. If this register is in the state 
|1), the previous position of the particle was n — 1. The 
coin register indicates if the walk should continue in the 
previously chosen direction (transmission in state |0)) or 
change the direction (reflection in state |1)). 

Taking into account the above, we define a shift oper- 
ator for a quantum walk with a 1-step memory on cycle 
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with d nodes as 

d-l 



Si = ^(|0)(0|<8)|0,n-l(modd))(0,i 



n=0 



+|0)(0|<8|l,n + l(modd))(l,n| (2) 
+|1)(1| (8 |l,n + l(modd))(0,n| 

+ |1)(1| ® |0,n- l(modd))(l,n| N 



or in a more consistent form as 
#1 = } " |c)(c| ® |/i m , c ,n + 2h miC - l(modd))(m, n|, 



E 

n.m.c 



(3) 

where /i m , c = m + c(mod2) represents a history depen- 
dence of the walk. 

The walk operator is defined as toss-a-coin and make- 
a-move combination, i.e. 

W 1 = S 1 {C®l 2 ®l d ), (4) 

where C is a coin matrix, e.g. Hadamard matrix H 

1/11 



yn 1 - 1 



(5) 



or any matrix C £ SU(2). 

The walk starts in some initial state \<fio). After each 
step the state is changed according to the formula 



\<t>n+l) = W?\<f>o) 

or as a recursive relation 

\<f>n+l) = Wi|^n). 



(6) 



(7) 



The probability of finding a particle at position v after 
n steps is obtained after averaging over the coin and the 
memory registers 



p(v,ri) 



c.m 



\(c, m,v\<l> n }\, 



(8) 



or, in other words, by tracing out over the memory and 
the coin subspaces 



p(v,n) = tr Cjm | (/>„)(</>„ |, 



(9) 



where tr c m denotes the operation of tracing out with 
respect to the coin and the memory subspaces. 

In |10) it was shown that p(v,n) is quasi-periodic for 
memoryless walks on cycles and it was suggested to con- 
sider quantity 



1 * 

s=l 



(10) 



which converges with t — > oo to the limiting distribution 
p(v). As the parameter n in this formula corresponds 
to the time required to perform n steps, we refer to so 
defined p(v) as time-averaged limiting distribution. 



III. PROBABILITY DISTRIBUTION 

In what follows we evaluate the probability distribu- 
tion of finding the particle at each node of the cycle. We 
consider the Hadamard walk only. Thus the walk opera- 
tor is given as 



W = S{H®t 2d ). 



(11) 



The second factor in the Eq. [TT] can be written in the 
matrix notation as 



H <g> l 2d = 



1 ( lid hd 

^2 \hd -hd 



(12) 



A. Fourier analysis 

In order to obtain expression for the probability distri- 
bution on the cycle we use Fourier analysis method (131 — 

Eg. 

Below we calculate amplitudes for a walk on cycle with 
a 1-step memory with the Hadamard matrix acting on 
the coin register. In this case we represent the vectors of 
amplitudes as 



$(n,t) 



/(0,0,rW 
(0,1, n|&) 

(1,0, n|&) 

\(l,l,n|&)y 



(13) 



The shift operator is defined as in Eq. ^ . In this case 
the interesting part of the shift operator reads 

S x (n) = |0)(0| ® |0,n)(0,n + l| + |0> <0| ® |l,n)(l,n- 1| 

+ |1)(1|<8> |0,n)(l,n + l| + |1)(1| ® |l,n)(0,n-l|' 

(14) 



After one step of time evolution we have 
$(n,t+l) = 



/{0,0,n+l\(H®l 2 ®ldMt^ 
(0, 1, n - l|(.ff ® 1 2 ® ld)\i>t) 
(l,l,n + lj(JT®l 2 ® ld)\ilt) 
\(l,0,n-l|(ff®l 2 ®ld)W) 



Evaluating the action of the Hadamard gate on the coin 
register one gets 



1 



' <0,0,Ti+l|i/j t ) + <l,0,n+l|Vt) ' 
$(»l,*+l)= -j= I ( ,l,„+l|V. t }-<l,l,n+l|Vt) 

(0,l,n-l|j/- t >+(l,l,n-l|V>t) 
(0,0,n-l|V>t>-(l,0,n-l|'/'*) I 



1 

71 



Rewriting the above expression using $(n + l,t) and 
$(n — 1, t) one gets 

$(n,t + l) =M-$(n + l,t) + M + ${n-l,t), (15) 
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where M + {advancing) and M_ [retarding) matrices read 

(16) 
(17) 



1 /0 s 

M — (0101 

M + 7= I | i 

\l I \ 1 -1 , 



1 / 1 1 
M - f 

7W - - I o i o -l 

\0 



In order to obtain the expression for the amplitudes 
of the quantum walk with memory on cycle we use the 
method introduced in |13j and represent time evolution 
of the walk using the Fourier transform 



d-l 



B. Initial state and time evolution 

In order to evaluate the probability distribution we 
need to choose the initial state of the walk. Again fol- 
lowing [T3] we start from the position with the coin in 
the state |0). As in this case we have extra register for 
storing memory, in the initial state the state vector reads 



(21) 



$(*, t + 1) = e 2vikn/d $(n, t + 1) (18) 

d-l 

= e 2wlkn/d (M_<f>(n + l,t) + M+$(n - l,t)) 

n=0 

From the above we get a recursive relation for the time 
evolution in the Fourier basis 



which means that the memory register is in the superpo- 
sition ^75 (|0) + |1)) and the coin register is in the state 
|0).The vector of amplitudes in this situation reads 



$(n,0) 



1 



0(n,t) = Afjt*(*,O), 



(19) 



where 



M, 



2ikn 
d 





e d 


\ 





2ifc7T 

e d 





2ifc7T 

e d 





e d 





2ifc7T 

— e d 


2ik-K 
e d 





—e d 


/ 



(20) 



V2\°oJ 

Hence there is only one term in its Fourier transform 



(22) 



(23) 



and thus in the Fourier basis we start from the state 



$(fc,0) 



1 (\ 



V2V0 



> 



(24) 



for any k = 0, . . . , d — 1. 

C. Initial state decomposition 



Let us now denote ^4 = exp(— 2irik/d). Matrix Mk has 
the following eigenvalues 

Ai = -1, 
A 2 = 1, 

1 + A 2 + VA 4 - 6A 2 + 1 



A, = 



2 V / 2A 



_ 1 + A 2 - \M 4 - 6A 2 + 1 
A4 " 272^1 ' 

with corresponding (orthonogal, but unnormalized) 
eigenvectors 



!'3 



2A+V2 

1 

V2A + 1 
A(^/2A + 2) 
2A + V2 
1 



> ) ^2 



V2A-1 
1 

A(v / 2A-2) 
V2-2A 
1 



Vi = 



±(A 2 -VA*-6A 2 + l-l) 

A 2 + V A 4 -6 A 2 + 1-1 
2A 2 



i(A 2 + V-4 4 -6A 2 + l-l) 

-A 2 + \/A 4 -6A 2 +l + l 
2A 2 
-1 
1 



Using the eigendecomposition of the matrix Mk we can 
calculate the form of the amplitudes in the Fourier basis 
after t steps. After t steps the vector of the amplitudes 
reads 



$(fc,t) = M^(fc,0) = M^ Q 



(25) 



for any k. 

Let us now write the initial state of the walk $(fc, 0) 
in the {«i}i=i,...,4 basis as 



where 



${k,0) = Y,a i (k)v i {k). 



( (fc) = (« i (fe),#(fe ) 0)) 



(26) 



(27) 



are components of $(fc,0) in {vi}i=i,...,4 basis. The evo- 
lution in the Fourier basis can be now written as 



$(M) = m£ £ <*(*)«<(*) = 5>i( fe )W«i(*0- 



(28) 
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The original components can be expressed by the 
Fourier-transformed vectors as 



d-l 



a — ^ 



fc=0 
-1 4 



(29) 



k=0j=l 



Using the above, the probability of finding the particle 
at n-th noted after t steps reads 

p(n,t) = Mn 7 t)\ 2 

d-l 4 

E E e 2ffi(m - fe)n/d <4(fc) ai (m) (30) 



xvt(fc)v,(m) [A J -(fc)*A i (m)] i 



D. Time-averaged limiting probability distribution 



Following [10 , let us now consider time-averaged prob- 
ability distribution p(n) = lim^oo \ X)l=oP( n ' s ) 



the sums are finite, we get 

d-l 4 



Using the expression (30) for p(n,t) and the fact that 



P( n )= E E K n (k,j,m,l)]im ~E[A}(*)AiM]' 

(31) 



fc,m— J,Z=1 



s=0 



where 



K n (k,j,m,l) = ^a*{k)a l {m)v]{k)v l {m)e 2 ^ m - k ^l d . 

(32) 

One can observe that the convergence of p(ri) depends 
only on the behavior of the term 

1 ^ 

f{k,j,m,l)=\\m - V[A*(fc)A ; (™)] s . (33) 

The value of this functions depends on the product of 
eigenvalues as 

f(k 7m/)-/ 1 if A K fc ) A '( w ) = 1 (34) 
I\K,3,m,i) < Q otherwige ■ W 



Definition in Eq. 34 can be easily calculated using the 
standard computer algebra systems, and thus allows for 
the evaluation of Eq. [31] 

The examples of time-averaged limiting distribution 



for the discussed model calculated using Eq. 31 are pre- 
sented in Fig. [T] In order to illustrate the influence of 
the initial state on the resulting distribution, two initial 
states were used to obtain these plots. In the case of the 
input state given by Eq. (21 1, the memory register is in 



the superposition, and the resulting time-averaged proba- 
bility distribution is symmetric with respect to the start- 
ing position. On the other hand, if the memory register is 
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FIG. 1. Limiting probability distributions for Hadamard 
quantum walks with one-step memory on cycle with d = 
3,4,6, 11 nodes for the initial states given by Eq. (21 1 (solid 
red line) and = (1,0,0,0) T (dashed black line). For each 
d the results are plotted in range [0, . . . , d] to illustrate the 
periodicity of the limiting distribution. Nodes are numbered 
starting from 0. 
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set to |0) in the initial state, the resulting time-averaged 
probability distribution does not have this property. 

The important difference in comparison with the mem- 
oryless quantum walk on cycle is that the time averaged 
limiting distribution depends on the initial state of the 
coin and the memory registers, but not on the parity 
of the number of nodes d [TTJ [TB]. The dependency is 
expressed in the otj(k) coefficients, which are defined in 
Eq. p7|. 



also pointed out the most important differences between 
quantum walks on cycles with and without memory. The 
most important difference is that the symmetry of the 
time-averaged limiting probability distribution is inde- 
pendent of the parity of the number of nodes. However, 
this distribution is heavily influenced by the initial state 
of the memory register. 
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